Possible three-dimensional chiral charge ordered superconducting state in cuprates 
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The 2D pair-condensate is characterized by a fluctuating chiral charge ordered state with a "checker- 
board" pattern in the CuO-z planes and with an alternating supermodulation along the c-axis in such 
a way that the adjacent layers are mirror images of one anothers electronic state. Planar chiral order 
is revealed with a recent circular dichroism (CD) ARPES experiment (A. Kaminski, et ai, Nature, 
416, 611. (2002)). We propose further CD experiments on ultrathin films with varying thickness 
and argue that the odd number of unit cells along the c-axis might provide dichroism hence might 
support the picture of 3D chirality in cuprates. We find that Coulomb energy gain occurs along the 
c-axis within a multilayer chiral charge ordered state, which is proportional to the measured bilayer 
condensation energy and to T c at optimal doping. Within our approach the superconducting (SC) 
pair is composed of the hole content of the coherence area and the self-repulsion of the condensate 
is compensated by the gain in the inter-layer Coulomb energy below T c . The SC condensate and 
the fluctuating charge order can also be described by a dynamical inter-layer electrostatic comple- 
mentarity. 
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I. INTRODUCTION 

A series of recent experiments explored the emergence 
of spin and charge order in cuprates [1-8]. Local probes 
such as scanning tunneling microscopy (STM) revealed 
recently evidence for a static checkerboard (CB) charge 
pattern with a real-space modulation periodicity of 4ao 
in the vortex core of Bi2212 [5], which is a provocative 
evidence for pinned charge stripes [6] . A 4ao X 4a charge 
pattern has also been reported in the absence of the ap- 
plied magnetic field [8], which is possibly induced by 
impurities at the surface. The STM observations could 
be consistent with a dynamical fluctuation of charge or- 
der states (COSs) which is expected to slow down with 
the increasing magnetic field and a static charge ordered 
pattern emerged in and around the vortex core [10]. It 
has been suggested that superconductivity and magnetic 
charge order coexist at moderate magnetic field and spin 
density wave order grows with the applied magnetic field 
[4]. Static magnetic charge order has also been observed 
in Y Ba2CuzOi-8 (YBCO) by neutron scattering exper- 
iment without the need of impurity doping and seems 
to be compatible with superconductivity as long as spin 
order remains dynamic [2] . The ground state of cuprate 
superconductors is also characterized by the interplay be- 
tween competing and coexisting ground states [7]. A 4- 
unit-cell superstructure (4ao x 4ao) is also reported in a 
high-energy X-ray diffraction study in optimally doped 
YBCO [9]. 



The CB COS of cuprates attracted recently the at- 
tention of several theoreticians as well [10-12]. It has 
been proposed that the magnetic field destroys the phase 
coherence of the hole pair by localizing them into a crys- 
tal, without breaking the pair [10]. In particular the CB 
modulated superconducting (SC) COS has been found fa- 
vorable in an intermediate doping interval using the t-J 
model [13]. It is still a question whether the modulations 
are present only in the local density of electronic states 
(LDOS) or also in the charge density [14,15]. If the SC 
COS can be described by a weak CB modulation [14], 
little modulations of the charge density may be hard to 
detect in the integrated LDOS [5] hence the improvement 
of the experimental resolution will be important in the 
future. 

Following a recent communication [16], in this work, 
we propose a specific microscopic state for the SC state 
in cuprates. This state has a 4ao X 4ao CB symmetry as 
observed in the experiment together with a local chiral 
symmetry (inter-layer mirror symmetry). A phase tran- 
sition to this kind of a state where two-dimensional (2D) 
parity (P) and time-reversal symmetries (T) are simul- 
taneously violated has been proposed and studied in a 
number of articles [17,18]. Up to date only planar chi- 
ral order is found [19] and early circular dichroism (CD) 
experiments have lead to controversial results in the ob- 
servation of three-dimensional (3D) chirality in cuprates 
[20,21]. We propose to recall these experiments and sug- 
gest further measurements which might support the ex- 



istence of a possibly hidden (local) electronic 3D chiral 
order in the superconducting state based on a number of 
a theoretical reasonings and experimental findings (sec- 
tion IV). 

We propose a simple phenomenological model for ex- 
plaining HTSC in cuprates in which these materials are 
nearly in a real space pairing limit [26]. The pairs of 
holes form two particle bound states, and then Bose con- 
dense at T c with the spatial constraint that the number of 
pairs per coherence area is approximately equal to 1. In 
order to avoid the frustration caused by the propagation 
of hole pairs [26] we propose a " scattered" density wave 
model of the hole pair state, in which the charge den- 
sity of the holes are distributed in the coherence area. 
In this sense the CB pattern might be an interference 
pattern due to elastic quasiparticle (QP) scattering on 
the Fermi surface of the superconductor [5] . Inter-layer 
(IL) Coulomb energy gain may arise due to the occu- 
rance of an alternating (asymmetric) interference charge 
pattern in the adjacent layers. In this case the sign of 
the broken P and T symmetry alternates between planes 
and the broken symmetry may be detected in principle 
by surface-sensitive probes or by weak effects in neutron 
scattering [27]. 

We would like to show that instead of using a com- 
plicated quantum formalism, it is possible to understand 
certain properties of HTSC dealing only with a classical 
picture of inter-layer electrostatics of cuprates. Although 
for simplicity we present a static charge order, the real 
charge order is dynamic in the SC state (with static spin 
density wave order). The static charge order might only 
appear upon high external magnetic field (together with 
dynamic spin order), which breaks, however down SC due 
to the pinning of holes to the lattice (COS of the Wigncr 
solid). Nevertheless, we model here the SC state by a 
static COS which can be taken as a "snapshot" of the 
dynamical picture. The resonance of these pictures leads 
to the realistic picture of the dynamical charge order in 
the SC state (charge fluctuations). 

Furthermore we assume the alternation (supermodu- 
lation) of the "checkerboard" charge pattern along the 
c-axis (that is normal to the planes) which leads to 
Coulomb energy gain (see FIG 1). We will discuss 
the possible experimental support of this assumption 
(see section IV). Without this assumption inter-layer 
Coulomb instability (CI) occurs in layered cuprates due to 
the enormous IL repulsion of holes. At best of our knowl- 
edge there is no widely accepted theory for avoiding the 
enormous IL Coulomb repulsion induced by the holes, al- 
though the role of holes in superconductivity is essential. 
In any naive model of electron pairing in cuprates the in- 
plane Coulomb repulsion is also troublesome. When the 
pairs of charge carriers are confined to the sheets, natu- 
rally a net self-repulsion of the pair condensate occurs. 
Although short range Coulomb screening in the dielec- 
tric crystal can reduce the magnitude of the repulsion 
it is insufficient to cancel Coulomb repulsion completely 
[30]. We would like to show that using a simple CB COS 




FIG. 1. The alternating chiral "checkerboard" charge or- 
dered state of the bilayer hole-antihole condensate mirror im- 
ages in the 4oo x 4ao bilayer superlattice model. Only a bilayer 
is plotted, however, this bilayer pattern alternates along the 
c-axis as well as the pattern is periodic along the ab-plane. 
Each lattice sites (opened (holes) and filled (antiholes) circles) 
correspond to a CuC>2 unit cell. Note the charge asymmetry 
between the adjacent layers. The bilayer can accommodate a 
pair of a boson condensate (4e). The inter-layer charge com- 
plementarity of the charge ordered state is crucial for getting 
inter-layer Coulomb energy gain. 

picture one can avoid both the in-plane and the out- of - 
plane CI of the layered SC condensate. 

II. CHARGE DISTRIBUTION IN THE 
CHECKERBOARD CHARGE ORDER STATE 

The CB COS can be understood as the in-plane al- 
ternation of hole-electron pairs in such a way that the 
Cooper wave- function is " composed" of the hole content 
p of the coherence area. The lattice site centered partial 
charges of the QP charge density is denoted as antiholc 
partial charges (q a h)- The sum of antihole charges at var- 
ious lattice sites correspond to the 2e charge of the quasi- 
particle (QP) within the coherence area (charge sum rule: 
J2i Qah,i ~ 2e, where q a u,i is the partial antihole charge 
at lattice site i). Therefore within our representation of 
HTSC the QP charge density is shared by each CuC>2 
sites in the coherence area. This is a marked difference 
with other theories (e.g. RVB [22]) where the SC pairs 
are formed between electrons of neighboring sites. In 
these theories the entire coherence area does not take 
part in the formation of QPs and hence the problem of 
self-Coulomb repulsion is huge. 

In the CB representation of the SC state every sec- 
ond lattice site {CuOi sites) is occupied by antiholes 
(the scattered nodes of the QP) and the rest are holes 
with partial charges of q^. 25 CuOi unit cells can be 
found in the 4a x 4a CB, and 25 x p ~ 4e hole con- 
tent is provided by the coherence area at optimal dop- 
ing, where po « 0.16e is the optimal hole concentration 
found in various single- and double-layer cuprates [23] 
which are described by the empirical relation T c /T™ ax = 
1 — 82.6(p — 0.I6) 2 [23] where p, the hole content varies 



in a broad range, pa « 0.16e is also found by first 
principles calculations as an optimal hole concentration 
[24]. Although p is measured in the normal state (NS) 
of cuprates, one can assume, based on simple chemical 
and physical intuity, that the chemical affinity of the 
CuC>2 unit might be similar in the SC and in the NS 
of cuprates. Therefore, if p 0.16e in the NS, similar 
quantity can also be assumed for the condensate anti- 
hole charge/CuC>2 units. For simplicity we consider only 
optimal doping, extension of this speculation is straight- 
forward for the entire doping regime by varying q a h,i ex- 
pecting a correspondence between p and the QP partial 
charges in the SC state. 

When the hole content condenses below T c , 2e charge 
is used for the charge neutralization of every second hole 
sites leading to an antiferromagnetic (AF) background 
COS possibly also with a CB modulation. The SC pair 
is formed from the other 2e charge coming from the 4e 
hole content on the AF background. The reason for the 
separation of the 4e hole content of the coherence area 
into 2e charge neutralizing charge and to 2e QP pair is 
a somewhat hypothetical, however seems to be a reason- 
able choice. By virtue only a single QP pair can be local- 
ized within the coherence area and therefore the other 2e 
charge is strongly coupled (pinned) to the lattice forming 
a delocalized AF background. 

The QP charge condenses to the charge neutralized 
sites leaving there antihole excess partial charges. The 
partly charge neutralized hole doped sheet possesses also 
a CB pattern with an antiferromagnetic (AF) spin or- 
der. The QP 2e charge condenses to the neutralized sites 
forming the antiholes in this AF background and a SC 
CB charge pattern occurs. The excess charge density 
at the antihole sites leads to the increase of the on-site 
repulsion. In order to maintain the phase coherence of 
the QP nodes (antiholes), the pairing-glue is necessary to 
compensate the effect of the increased on-site repulsion. 
We argue in the rest of the article that the pairing-glue is 
provided by the gain in the IL Coulomb energy and which 
supports pairing up to T c where the thermal fluctuation 
breaks the phase coherence of the QP pair. 

Charge fluctuations (CF) might occur due to the hop- 
ping of the antiholes to the holes (hole annihilation) . CF 
is a collective phenomenon in the AF background (the 
coherent motion of holes), the antiholes hop from an an- 
tihole site to a hole simultaneously keeping the phase 
coherence of the 2e pair and leading to a fluctuating or- 
der [6]. The charge fluctuations are associated with an 
electronic liquid crystalline state (melted Wigner crys- 
tal) [14]. This kind of a dynamic picture of the charge 
density wave order might be consistent with the presence 
of a periodic modulation in the electron hopping or pair- 
ing amplitude [15]. Therefore the collective hole-antihole 
"exchange" between neighboring lattice sites leads to a 
resonance between CB COSs. Within our picture the 
motion of a hole is accomplished by the coherent motion 
of partial holes (the nodes of the scattered hole pair) as 
well as the motion of a QP is associated with the co- 



herent motion of partial charges (antiholes, the scattered 
nodes of the condensate charge density waves) in the liq- 
uid condensate. The site-by-site motion of a complete 
hole (a complete missing electron at a site) is an unfa- 
vorable process with a decent barrier and leads to an 
energetically frustrated system [26]. Instead we propose 
the motion of the scattered density waves or holes on the 
lattice. 

Distributing the 2e charge of the quasiparticle conden- 
sate in the 4oq x 4ao coherence area the black and white 
" fields" of the CB correspond to lattice site centered an- 
tiholes and holes with partial charges of qh,ah ~ 4e/A 2 w 
±0.16e, where N is the real space periodicity of the co- 
herence area (N « Cab/ a o + 1,N = 5 for the 4a x 4a 
lattice, iV 2 is the number of CuOy, unit cells in the co- 
herence area) where £ a {, « 16 A is the in-plane coherence 
length [25] and a w 3.9A is the lattice constant. 

We would like to calculate the IL electrostatic energy 
of an alternating multilayered CB COS, which is essen- 
tial to avoid Coulomb instability even in the vortex core. 
The multilayered CB COS might be energetically com- 
petitive with other possible non-alternating pattern in 
the SC state (e.g. stripes). In the CB state the conden- 
sate does not " suffer" from the in-plane and out-of-planc 
CI, and also energy gain occurs due to the out-of-planc 
alternation of holes. Within the planes the net Coulomb 
self-interaction nearly vanishes due to the equal number 
of repulsive and attractive pair interactions (inter-site) in 
the CB alternation of holes and antiholes. In other words, 
the planes are sufficiently relaxed from Coulomb forces 
(electrostatically balanced), otherwise the 2D structure 
could not be retained. Indeed, in most of the cuprates 
no out-of-plane distortion of the CuOi layers (buckling) 
is observed, except in YBCO. Even in this case buckling 
is weak and is observed in the NS [33]. Other possible 
COSs (e.g. stripes) with no CB pattern have a tendency 
to be electrostatically unbalanced and, therefore, ener- 
getically unfavorable. The only source of the energy gain 
in the potential energy is provided by the IL Coulomb 
coupling which is converted to the condensation energy 
of the SC state (the free energy gain when the system 
goes to the SC state). 

The kinetic energy driven mechanism of HTSC is also 
proposed as a possible source of pairing [26,28]. It must 
be admitted that we can not rule out the kinetic energy 
gain mechanism and more sophisticated approaches will 
be necessary to get a clear cut conclusion on this prob- 
lem. We simply remark here that owing to the law of the 
conservation of energy the gain in the potential energy 
side must be converted to an increase in energy on the 
kinetic side and vice versa. This basic law implies that 
even if saving in kinetic energy is responsible for HTSC 
then there must be an energy loss on the potential energy 
side. The advocates of the kinetic energy mechanism at- 
tribute the potential energy loss naturally to the pairing 
induced increase in the planar repulsion energy [26]. In 
this paper we would like to show that it is possible to ac- 
count for T c and for the condensation energy using only 
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FIG. 2. The 3D view of the static "checkerboard" charge 
pattern of the boson condensate corresponding to the 
4ao x 4ao charge modulation represented by the order pa- 
rameter given in Eq. (1). x and y coordinates are given in 
A. The wells and peaks correspond to holes and antiholes, 
respectively. 



the charge complementarity induced IL Coulomb energy 
gain. 

We would like to study then the magnitude of direct 
Coulomb interaction between charge ordered square su- 
perlattice layers as a possible source of pairing interac- 
tion. Our intention is to understand HTSC within the 
context of an IL Coulomb-mediated mechanism. The IL 
charging energy we wish to calculate depends on the IL 
spacing, the IL dielectric constant e c , the hole content p 
and the size of the superlattice (£ b)- Finally we calculate 
the static c-axis dielectric constant e c and the coherence 
length £ a b for various cuprates which are compared with 
the experimental observations. 

III. THE SUPERCONDUCTING ORDER 
PARAMETER 

The superconducting order parameter (OP) of the d- 
wave condensate which corresponds to the model with a 
checkerboard charge modulation in the planes takes the 
form of 

y) = nl /2 [cos(— ir) + cos(—ir)]. (1) 
ao ao 

For simplicity the distribution of the OP is neglected 
in the 3rd dimension and a nearly perfect 2D character 
is attributed to the condensate. The 3D anisotropy of 
the condensate is negligible in the superconducting (SC) 
state which is reflected by the ratio of the in-plane and 
out-of-plane coherence lengths £, a b/£,c ~ 10 [25]. The fac- 
tor no is the maximal value of the charge density at the 
lattice site centers. Eq. (1) is displayed for the coherence 
area in FIG 2. The modulation of the order parameter 



corresponds to the real-space modulation of the hole den- 
sity in the superconducting (SC) state. This kind of an 
order parameter is given earlier by Alexandrov [11]. 

The order parameter must satisfy the charge sum rule 
for the boson condensate indicating the localization of 
the pair condensate within the coherence area, 

rial, 

2« / \*{x,y)\ 2 dxdy. (2) 
Jo 

Another restriction on ^(x,y) is that its integral over a 
unit cell with the area of ~ (ao/2) 2 must correspond to 
the partial point charges of holes and antiholes in the CB 
COS, 

/■ao/2 

\lh(ah)\ = / \*{x,y)\ 2 dxdy. (3) 
Jo 

Qh(ah) can a l so be taken as a free parameter ("filling" fac- 
tor) at the lattice sites. We assume homogeneous filling, 
all the lattice sites occupied by the same q a h or Cal- 
culations indicate that the magnitude of the calculated 
condensation energy is in accordance with the experi- 
ment when quah) ~ P ~ 0.16e. The variation of qh( a h) 
introduces the effect of doping. This equation reflects 
the lattice site centered localization of the holes and an- 
tiholes and leads to a simple electrostatic model where 
the charge modulation for simplicity is replaced by clas- 
sical point charges centered in the center of the holes and 
antiholes. However, it must be admitted that in certain 
cuprates, usually with low T c , the coherence area is larger 
then 4a x 4a . In these cases qu(ah) ^ 0.16e at opti- 
mal doping. For instance in LSCO (Lai ^Sro.i^CuO^s) 
S^ab ~ 7ao, therefore qh( a h) — 4e/64 w 0.063e at opti- 
mal doping which leads to a weak modulation both in 
the LDOS and in the charge density representing a great 
challenge for the experiment. 

One can also speculate on the exponential decay of the 
OP at the border of the coherence area leading to inho- 
mogeneous distribution of quah) within the CB similar 
to that of in refs. [10,12]. For simplicity we are dealing 
only with the homogeneous distribution of qh(ah) within 
the CB which can be taken as an average distribution of 
the charge pattern. Charge inhomogcncity could also be 
troublesome since partial charges \qh( a h) \ > 0.16 could be 
energetically unfavorable due to the increasing on-site re- 
pulsion with increasing charge density at certain regions 
of the CB. 

IL Coulomb energy gain occurs only in that case when 
holes in one of the layers are in proximity with antiholes 
in the other layer (FIG 1 , IL electrostatic complementar- 
ity, bilayer model). An important feature is then that 
the boson condensate can be described by an IL charge 
asymmetry (mirror symmetry). Therefore we assume an 
alternating charge pattern along the c-axis. The IL cou- 
pling of the boson-boson pairs in the bilayer 5x5 model 
naturally suggests the effective mass of charge carriers 
to* w 4m e , as it was found by measurements [29]. Fur- 
thermore, the CFs in different planes are also coupled, 



in order to keep charge complementarity even during the 
fluctuation (" superexchange" ) of the CB charge patterns. 

Important to note that the charge sume rule holds 
for the characteristic bilayer with a coherence area 

1i hole — 4e. In other words a pair of a boson con- 
densate can be localized within a characteristic bilayer 
depicted in FIG 1 . The alternation of the charge pattern 
along the c-axis is reflected in the order parameter by 
the layer-by-layer alternation of cos and sin functions as 
well as the in-plane CF might also be reflected by the 
" superexchange" (coherent self-exchange) of cos and sin 
functions due to the collective motion of holes and an- 
tiholes. This is nothing else then the manifestation of 
the sign reversal of the order parameter which is seen by 
tunneling experiment as well [31]. 

IV. POSSIBLE EXPERIMENTAL 
FINGERPRINTS OF THE CHIRAL ORDER 

Although a direct experimental evidence is not re- 
ported until now for the c-axis alternation of the hole 
density we discuss certain measurements which might be 
related to the presence of a chiral order in cuprates. The 
experimental detection of the electronic chiral order could 
be difficult in the geometrically achiral cuprates (due to 
the lack of the geometric inversion center, noncentrosym- 
metric compounds) therefore we discuss various measure- 
ments to get a comprehensive picture. 

The recently published circular dichroism ARPES ex- 
periment [19] is consistent directly with a planar chi- 
ral order. The experimental setup in these experiments 
incorporates a mirror plane perpendicular to the de- 
planes hence the chirality of the ab-plane is explored 
((0, 0) — (tt, ir)(T — M line)) [34]. This measurement pro- 
vides a provocative evidence for a broken symmetry state 
below T*, where a pseudogap appears in the underdoped 
Bi2212 compound. Originally it was proposed by Varma 
et al. to support the circular orbital current model [35] 
which might be consistent with a circular charge mo- 
tion picture within a CB pattern with a coherence area. 
According to our proposal with IL coherence of charge 
motion the circular currents are must be in the same 
direction with each other in the adjacent layers to ob- 
serve dichroism also along the c-axis. Circular currents 
with opposite orientation in the adjacent layers would 
not lead to dichroism because the adjacent layers form 
an enantiomorphous pair with globally restored symme- 
try hence no difference in absorption between right and 
left circularly polarized radiations could be found. 

Unfortunately an experimental test of 3D chirality is 
difficult to make directly. It is hard if impossible to re- 
peat the CD ARPES measurements with a mirror plane 
parallel to the a6-plane owing to the fact that thin films 
are in general use in ARPES measurements [34] . Another 
problem is that in many AF crystals, where the locally 
time-reversal symmetry can be globally restored, the to- 



tal dichroic signal is zero [36]. The failure of observing 
circular dichroism in cuprates [20] could be attributed 
then to the macroscopic symmetry of these systems (hid- 
den chiral order), for example to the formation of an op- 
tically inactive phase (a racemate with enantiomorphous 
pairs of layers). Globally no dichroism occurs also in 
the c-axis alternating CB COS superstructure because 
the crystal is stacked by enantiomorphous pairs of lay- 
ers (also without the assumption of circulating currents). 
We place in this case the inversion center to the midpoint 
of the coherence bilayer shown in FIG 1 . 

Ultrathin films with odd number of layers, however, 
should provide nonvanishing dichroic signal (optical ro- 
tation) because of the presence of an unpaired optically 
active layer. In ultrathin layers the number of layers can 
be kept under control, even one-unit-cell thick sample can 
be prepared [40] . Its signal is not compensated by a signal 
of its mirror image hence should be detected by CD mea- 
surements. Therefore, we propose to repeat the early CD 
measurements with odd number of layers, hence double- 
layer cuprates (such as YBCO or Bi2212) are unsuitable 
for the detection of 3D chirality with traditional CD tech- 
niques. The CD measurements on thin layer YBCO [20] 
therefore naturally does not result in 3D chirality. In- 
stead the ultrathin film of the single-layer LSCO should 
be used with odd number of unit cells along the c-axis. 
Although Lyons et al. [21] found no sign of optical activ- 
ity in LSCO, the application of specific conditions pro- 
posed above might indicate optical activity in the under- 
or optimally doped LSCO. Ultrathin films with varying 
thickness made from the three-layer BH223 compound 
should also be suitable for CD studies. Within our model 
a one-unit-cell thick single-layers do not show HTSC due 
to the lack of IL coupling. Gelfand and Halperin [39] 
have also proposed the refinement of optical experiments 
with odd number of layers which would reveal clear ev- 
idence of broken time-reversal symmetry and of anyons 
in cuprates. At best of our knowledge since then no ex- 
periment has been done using such conditions. 

Moreover, the dichroism signal should alternate with 
the parity of the number of layers in the thin films: the 
even number of layers should give vanishing signal while 
the odd number case should result in a decent absorption 
difference between the right and left circularly polarized 
radiations. It should be noted that we do not question 
the remark of Lawrence et al. [20] that the strong temper- 
ature dependence of the dichroic signal found in YBCO 
and in bismuthates [21] must be attributed to surface 
inhomogeneities and to contaminations. We simply pro- 
pose to apply these experiments for other compounds 
with the systematical variation of the film thickness to- 
gether with the careful treatment of the surface. The 
strong CD signal below T c found in the single-layer bis- 
muthate (Ki- x Rb x BiOz) [21] according to our knowl- 
edge is never falsified. 

In any case we also propose to probe the chiral order 
in the slightly underdoped cuprates in the presence of a 
strong magnetic field where a quasi static charge order 



emerged possibly with an internal asymmetry. In the 
SC state, however, the fluctuating magnetic order might 
weaken the signals coming from the presence of the chiral 
order hence the detection of the charge asymmetry could 
be difficult. 

Recently the anomalous X-ray diffraction has been pro- 
posed to circumvent the limitation of the global symme- 
try, where the local transition amplitudes are added with 
a phase factor that can compensate the vanishing effect 
due to the global symmetry [36]. Unfortunately no ap- 
plication of X-ray dichroism is known for cuprates until 
now. Optical activity has been detected using the natu- 
ral circular dichroism in quite a few noncentrosymmetric 
crystals [37] . The x-ray magnetochiral dichroism has also 
been proposed to unravel hidden space-time symmetry 
properties of magnetoelectric crystals [37] . The applica- 
tion of these techniques for cuprates might also reveal a 
possible hidden chiral order in the SC state. At present 
the experimental test of a chiral COS in SC cuprates re- 
mains, however, indirect. 

There are several indirect experimental findings which 
support the IL scenario, such as the T c variation in thin 
film heterostructures [40,41]. In such films HTSC van- 
ishes when the film thickness approaches the IL distance. 
The variation of T c upon the external hydrostatic pres- 
sure in Hg-cuprates compounds [32] can also be explained 
by the IL coupling concept, since upon pressure the IL 
spacing of cuprates decreases seriously which leads to an 
increase in T c due to the possible enhancement of IL cou- 
pling. 

The ARPES circular dichroism experiment reveals that 
the breaking of the time-reversal symmetry appears in 
the underdoped samples while is not present in the over- 
doped regime [19]. The absence of planar chirality in the 
overdoped samples might be consistent with a 2D <-> 3D 
charge transfer of the hole content [16] which does not 
occur in the underdoped regime when the planes remain 
coupled hence mirror symmetry is retained. In the over- 
doped side, however, the strong 2D <-> 3D charge fluc- 
tuations screen IL coupling hence IL symmetry progres- 
sively develops with overdoping. Therefore overdoping 
might destroy the chiral order (decouples the layers) and 
favors the occurrence of a symmetric Fermi-liquid phase. 
Dispersion along the T — Z line in ARPES spectra would 
be the indicative of 2D <-> 3D CF although it is hard 
to collect such a data from presently available ARPES 
results [42]. 

The observed supression of the pseudogap in the over- 
doped side of the SC dome might also be associated with 
the enhancement of 2D <-> 3D charge fluctuations. The 
pseudogap vanishes above the critical hole concentration 
p « 0.19 [44] which might be associated with the decou- 
pling of layers. The superconducting peak ratio (SPR) 
seen in the ARPES spectra of Bi2212 peaks at p w 0.19 
and drops for p > 0.19 [45] which also support that IL 
decoupling effects become dominant in the strongly over- 
doped regime. This is because SPR can be related (not 
directly) to the superfluid density hence the drop of the 



SPR signal intensity reflects the reduction of the car- 
rier concentration in the planes coexisting with the high 
doping level and is associated with the 2D <-> 3D charge 
fluctuations together with the lost of the 2D order in the 
strongly overdoped regime. First principles calculations 
indicate the saturation of the in-plane hole concentra- 
tion around p « 0.19 and support the presence of the 2D 
<-> 3D charge fluctuations [24]. These calculations reveal 
that the hole concentration in the planes reaches its max- 
imum when the dopant band is completely occupied, e.g., 
for a closed O shell [24]. This occurs at around p ~ 0.19 
and a strong c-axis anisotropy of the hole content devel- 
ops as demonstrated in the calculated LDOS [24] which 
also peaks at p « 0.19 [24] in accordance with the specific 
heat anomaly [38,44]. Interestingly the pseudogap line 
(T*) crosses the SC dome atpw 0.19 [44]. For p > 0.19 
(strongly overdoped regime) T* < ksTc and the hole 
content can freely 2D <-> 3D fluctuate. Thus we predict 
the occurrence of a 3D dispersion in the ARPES spectra 
along the out-of-plane direction r — Z in the Brillouin 
zone in the strongly overdoped regime and above T c . The 
strongest evidence for a coherent 3D transport is found 
in the overdoped T12201 [43]. The reported polar angu- 
lar magnetoresistance oscillations in high magnetic field 
firmly establishes the existence of a coherent 3D Fermi 
surface. This is consistent with our 2D <-> 3D fluctua- 
tion model of the hole content in the heavily overdoped 
regime and in the normal state. 

The issue of bilayer splitting seen in ARPES spec- 
tra [38,42] could also provide important ingredients in 
searching for the signatures of chiral order and inter-layer 
coupling. The two different component in the spectra is 
attributed to bonding and antibonding bands with oppo- 
site symmetry along the c-axis with respect to the mid- 
point between the two CuOi layers in the bilayer block 
of Bi2212 [42]. The overdoping induced appearance of 
a band with a strong c-axis anisotropy (due to 2D <-> 
3D CF) could explain bilayer splitting. This band might 
also provide a peak at (n, 0) and could be the signature 
of the broken chirality (the recovery of symmetry) in the 
normal state of cuprates. This issue is still highly con- 
troversial and further specific studies are needed which 
provide further informations on the 3D nature of HTSC. 



V. THE CONDENSATION ENERGY AND T c 

We calculate the condensation energy/ C11O2 (Uo) in 
order to show that at a reasonable choice of the IL di- 
electric constant e c the calculated Uo is typically in the 
range of the measured values. We will show in this sec- 
tion that this correlation between the measured and the 
calculated values manifest in a number of ways. 

First the condensation energy of a bilayer with a co- 
herence area in the planes (Uob) is given as follows 



f7 b = 2(n+l) 



a 



+ 1 



IL,SC 



(4) 




20 40 60 80 
critical temperature (K) 
FIG. 3. The bilayer condensation energy (Uob, meV) as a 
function of the critical temperature (K) at optimal doping. 
The straight line is a linear fit to the data. The slope of 
the linear fit is Uob/T c m fcs which is a strong evidence of 
Eq. (7). The error bars denote standard deviations estimated 
from various measurements of the condensation energy. 



where E^ L ' SC is the Coulomb energy gain in the SC state. 
Uq is the experimental condensation energy given per 
unit cell. Eq. (4) is generalized for multilayer cuprates 
introducing n. For single layer cuprates n = 0, for bilay- 

1 2 



ers n = 1, etc. The factor 



Lk + i 

a 



is the number of 

CuOi lattice sites in the planes (« TV 2 , within the co- 
herence area). Factor 2 is applied in Eq. (4) because we 
calculate the condensation energy of a bilayer. The IL 
Coulomb energy is in the SC state 



E, 



IL,SC 



e 2 Q 
47re e c 



where 



Q 



Ni N 2 (n) (k) 

EE 



m— 2 ij 



(k,l) 



(5) 
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where r^'^ is the inter-point charge distance and r tJ 



(fe,0 



> 



do,, where d]L is the IL distance (Cu0 2 plane to plane, 
i j). k,l represent layer indexes (k =/= I), and cfp 
are the point charges of holes and antiholes centered at 
CuOi lattice sites in the fcth and Ith layers. First the 
summation goes within the bilayer up to iV 2 then the IL 
Coulomb interaction of the basal bilayer are calculated 
with other layers along the c-axis in both direction (k = 
1,2). Ni is the number of layers along the c-axis. When 
Ni -> oo, bulk E* L < SC is calculated. 

The plot of U b (the bilayer condensation energy) 
against T c is shown in FIG 3 using only experimental 
data. Remarkably the data points of various cuprates 
with a variety of critical temperature fit to a line and its 
slope Uob/T c is the Boltzmann constant ks- The aver- 
age value we get is ks ~ 1.3 ± 0.2 x 10~ 23 J/K which is 
remarkably close to the value of ks = 1.38 x 10~ 23 J/K. 
In the rest of the paper we will present further evidences 



in order to show that the agreement might not be acci- 
dental. According to the correlation found between Uob 
and T c the following formula can be given using Eq. (4) , 



U ob = 2(n+l) 



Cafe 

a 



1 



U Q rj k B T c . 



(7) 



Therefore the bilayer condensation energy Uob can di- 
rectly be related to the thermal motion at T c . It must be 
emphasized that Eq. (7) is coming directly from the em- 
pirical relation found in FIG 3 and its theoretical deriva- 
tion is currently not available yet. This relation at a first 
look is certainly unusual and is completely different from 
the formula known in conventional superconductors [25]. 

We understand this correlation as follows: the SC pair- 
condensate is stable against thermal fluctuations up to 
T c , since the IL coupling energy E^, L ' SC k B T c . Hence 
above T c the phase coherence of the QP pair is broken 
and SC is suppressed. For T > T c , in the underdoped side 
pinning of the holes becomes stable against CFs, whereas 
in the overdoped side the hole content exhibits a strong 
3D anisotropy hence the SC phase vanishes. As far as the 
microscopic mechanism of pairing is concerned, we fur- 
ther argue that the IL Coulomb energy gain is converted 
to the stabilization of the fluctuating CB COS, which 
is although electrostatically is balanced (free of frustra- 
tion), however, the increased on-site repulsion at the an- 
tihole sites (due to the increased density, q a h ~ —0.16) 
is energetically unfavorable and is compensated by the 
IL energy gain up to T c . Above T c therefore the on-site 
repulsion of the condensate is no longer compensated by 
the gain in the IL Coulomb energy. The on-site repulsions 
of the QP charge density sum up to the self-repulsion of 
the scattered QP pair which is compensated then by the 
out-of-plane Coulomb energy gain (the pairing glue) be- 
low T c . The magnitude of the self-repulsion of the con- 
densate with CB COS could easily be checked by sophis- 
ticated ab initio or Hubbard-type approaches. Our rough 
estimate based on simple in-plane Coulomb energy cal- 
culations provides values in the range of ~ fcflT c . If the 
on-site repulsion of the condensate is compensated by the 
gain in the IL Coulomb energy then it follows that the 
optimal hole content po « 0.16 is more or less general for 
various cuprates with different T c just because the coher- 
ence area increases with decreasing T c hence the on-site 
repulsion is also decreasing with decreasing antiholc par- 
tial charges (Table I). Larger then 4ao x 4ao coherence 
area naturally suggests that q a h < po, therefore the re- 
maining part of the hole content does not contribute to 
the QP charge density and is rather pined to the lattice 
below T c . Above T c the phase coherence of the fluctuat- 
ing CB COS is lost via hole annihilation (slowed down 
CFs in the underdoped regime) leading to a frustrated 
AF spin order (spin glass) or via the 2D <-> 3D phase 
transition of the hole content. 

It is still also a possibility that only the phase coher- 
ence is broken and the in-plane fluctuation of the holes 
persists (at least partly) above T c leading to the pseudo- 



TABLE I. The calculated coherence length of the pair con- 
densate given in a using the experimental condensation en- 
ergies of various cuprates and Eq. (8) at optimal doping. 

T c (K) k B T c (mcV) Uo (^eV/u.c.) O^M C7M 



Bi2201 20 

LSCO 39 

T12201 85 

Hgl201 95 

YBCO 92 

Bi2212 89 

Hgl223 135 



1.6 

2.5 
7 

7.8 
7.5 
7.3 
11 



10" 

21 6 
100 ± 20 d 
80 - 107 e 

110 s 

95 s 

114 



~ 7 
~ 5 
~ 5 
~ 3 
3-4 

4 



5-8 c 

5' 

3- 4 h 

4- 6 1 
4 f 



dQ 3.88A, Uo is the measured condensation energy of various cuprates 
in fi cV per unit cell at optimal doping. a from [58], b U ss 2 J/mol 
from [50,51], c from [25], d [59], " Uo ss 12 - 16 mj/g from [52,60] 
and £ ab from [53], d Uo ~ 11 J/mol from [54], ' Uo f» 10 J/mol 
from [54] , ^ from [53] , 5 from [54] , ^ from [25] , 1 from recent mea- 
surements of Wang et al. , £ a b ~ 23A(~ 5 — 6ao ) [46] ,from STM images 
of ref. [5] £ a b ~ 4ao, ^^ Ic is calculated according to Eq. (8) and is also 
given in Table I and ls t ^ le measured in-plane coherence length 

given in ao ~ 3.9A. The notations arc as follows for the compounds: 
Bi2201 is Bi 2 Sr 2 Cu0 6 +s, LSCO {L ai . sb Sr . lb CuO i+ s), T12201 
(Tl 2 Ba 2 C'u0 6 ), Hgl201 (HgBa 2 Cu0 4+ s), YBCO (YBa 2 CuO r ) and 
Bi2212 is Bi 2 Sr 2 CaCu 2 O s+s . Hgl223 is HgBaCa 2 Cu 3 12+s . 



gap phenomenon. The incoherent motion of the charge 
carriers might also be then associated with the NS gap 
either within the planes and out-of-plane. This issue cer- 
tainly awaits further theoretical investigations. 

In order to test the validity of the empirical Eq. (7) 
we estimate the coherence length of the pair condensate 
derived from Eq. (7) and using only experimental data, 



£ab 



a 



2(n+l)U 



(8) 



The results are given in Table I as (in a unit) and 
compared with the available measured Cab P ■ The agree- 
ment is excellent which strongly suggests that Eq. (8) 
should also work for other cuprates. We predict for the 
multilayer Hg-cuprate Hgl223 (T c = 135 K) the conden- 
sation energy 114 /xeV/u.c. and £ a b w 4a which is in 
agreement with the measured value [53]) using the di- 
electric constant of e c « 61 for the Ca-spaced trilayer 
block. 

The expression Eq. (7) leads to the very simple formula 
for the critical temperature using Eq. (4) and a simple 
Coulomb expression for the IL coupling energy E I C L,SC 



T C {N, d, e c 



e 2 Q 



47re e c fcB 



(9) 



When Ni — > oo, bulk T c is calculated. T c can also be 
calculated for thin films when iVj is finite. e c can also be 
derived 



e 2 Q 



TABLE II. The calculated dielectric constant e c using 
Eq. (10) in various cuprates at the calculated coherence length 



£ a b of the charge 


ordered state g 


;iven in 


Table I. 






d(A) 


Tc(K) a 


>») 


e c 


exp 


Bi2201 


12.2 


20 


8 


9.9 


12 a 


LSCO 


6.65 


39 


7 


11.3 


23 ±3, 13.5 6 


Hgl201 


9.5 


95 


4 


6.5 




T12201 


11.6 


85 


5 


13.0 


11.3 C 


YBCO 


8.5 


93 


3 


19.4 


23. 6 d 



where ^^ lc is the estimated in-plane coherence length given in ao ^ 
3.9A unit, d is the Cu0 2 plane to plane intcr-laycr distance in A, T c is 
the cxperime ntal critical temperature. e c is from Eq. (10). e e c xp are the 
measured values obtained from the following references: a [55], b [57], 
or from reflectivity measurements [59], u> p « 55cm -1 [61], A c ~ 3^im 
[60], c from [59], d from reflectivity measurements: lj p Ri 60cm -1 [61], 
A c ps 0.9^m [48]. 



where a c-axis average of e c is computed when Ni — > oo. 

The calculation of the c-axis dielectric constants e c 
might provide further evidences for Eq. (7) when com- 
pared with the measured values [47,57]. In Tabic II 
we have calculated the static dielectric function e c using 
Eq. (10) and compared with the experimental impedance 
measurements [56,57]. e c can also be extracted from 
the c-axis optical measurements using the relation [59] 
e c (uj) = e c (oo)—c 2 /(ujpX 2 ), where eoo andw p are the high- 
frequency dielectric constant and the plasma frequency, 
respectively [62]. c and A c arc the speed of light and 
the c-axis penetration depth. At zero crossing e c (uj) = 

1 /2 

and Up = c/(A c e c (oo)). Using this relation we predict 
for the single layer Hgl201 the low plasma frequency of 



i 8 cm using e c = 
8 lira [60] using £ ab 



e c (oo) w 6.5 (Table II) and 
5ao [53]. 



4:TTe k B T c 



(10) 



VI. CONCLUSION 

Our primary result is that the boson condensate in 
the superconducting state can be described by an alter- 
nating " checkerboard" type of charge pattern along the 
c-axis which leads to inter-layer charge complementarity 
(chiral charge ordered state) and Coulomb energy gain. 
Our proposal is that this gain is converted to the conden- 
sation energy of the superconducting state, although the 
detailed mechanism of this process still remains unclear. 
Within our model the pairing glue is also provided by 
inter-layer coupling. 

The superconducting state can also be described by a 
fluctuating order hence the phase coherence of the quasi- 
particles can only be retained if inter-layer charge and 
phase complementarity emerged. This physical picture 
naturally explains the variation of T c system-by-system 
in various conditions (external and chemical pressure, 
multilayers, heterostructures etc.). Without the assump- 
tion of inter-layer charge complementarity the supercon- 
ducting state would "suffer" from an enormous inter- 



layer Coulomb instability (repulsion) which is certainly 
not the case. 

The charge chirality model explains the recent circular 
dichroism ARPES experiment in which the underdoped 
samples show dichroism below the characteristic temper- 
ature T* of the pseudogap. We believe that the modula- 
tion of the charge density along the c-axis theoretically 
provides the possibility of better understanding HTSC. 

Finally we attribute the breakdown of SC on the over- 
doped side of the superconducting dome to the destroy 
of the 2D order (the 2D — > 3D transition of the hole- 
content in the normal state) and to the emergence of an 
antiferromagnetic order on the underdoped side. We sup- 
port a picture which incorporates a fluctuating order in 
the optimally doped regime separated from the Fcrmi- 
liquid-like overdoped regime by a quantum critical point 
and by the suppression of charge fluctuations from the 
underdoped region. 
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